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Abstract 



We study the dimensionally regularized fermion propagator Dyson- 
Schwinger equation in quenched nonperturbative QED in an arbitrary co- 
variant gauge. The nonperturbative fermion propagator is solved in D = 
4 — 2e < 4 dimensional Euclidean space for a large number of values of e. 
Results for D = 4 are then obtained by extrapolation to e — > 0. The non- 
perturbative renormalization is performed numerically, yielding finite results 
for all renormalized quantities. This demonstrates, apparently for the first 
time, that it is possible to successfully implement nonperturbative renormal- 
ization of Dyson- Schwinger equations within a gauge invariant regularization 
scheme such as dimensional regularization. Here we present results using the 
Curtis- Penningon fermion-photon proper vertex for two values of the coupling, 
namely a = 0.6 and a = 1.5 and compare these to previous studies employing 
a modified ultraviolet cut-off regularization. The results using the two differ- 
ent regularizations are found to agree to within the numerical precision of the 
present calculations. 



I. INTRODUCTION 



Strong coupling quantum electrodynamics (QED) has been extensively studied through 
the use of the Dyson-Schwinger equations (DSE) |]-[|. In such studies it is unavoidable that 
the infinite set of coupled integral equations be truncated to those involving Green functions 
with relatively few external legs. This truncation has the consequence that one or more of 
the Green functions appearing in the remaining equations are no longer determined self- 
consistently by the DSE's and so must be constrained, for example, by known symmetries 
(including Ward-Takahashi identities (WTI) Q), the absence of artificial kinematic singular- 
ities, the requirements of multiplicative renormalizability (MR), and must be in agreement 
with perturbation theory in the weak coupling limit. Furthermore the gauge dependence 
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of the resulting fermion propagator should eventually be ensured to be consistent with the 
Landau-Khalatnikov transformation ||. 

What makes QED a particularly attractive theory to study with DSE techniques is 
that the coupled integral equations determining the photon and fermion propagators are 
completely closed once the photon-fermion proper vertex is specified. Up to transverse 
parts this proper vertex is in turn determined in terms of the fermion propagator by the 
corresponding WTI. Thus, the state of the art for this type of calculation consists of imposing 
the greatest possible set of constraints and constructing the most reasonable possible Ansatz 
for the transverse part of the vertex. While such a DSE approach to nonperturbative QED 
can never be an entirely first-principles approach, such as a lattice gauge theory treatment 
f|, it has the advantages that it is at times possible to obtain some analytical insights, there 
is no limit to the momentum range that can be studied and one is able to compare different 
regularization schemes. This latter property will be made use of in this paper. 

A number of discussions of the choice of the transverse part of the proper vertex can be 
found in the literature, e.g., Refs. 0-15J- We will concentrate here on the Curtis- Pennington 



(CP) vertex [|g|- |IT]JT6|j , which satisfies both the WTI and the the constraints of multiplicative 
renormalizability. With a bare vertex (which breaks both gauge invariance and MR) the 
critical coupling of quenched QED differs by approximately 50% when calculated in the 
Feynman and Landau gauge. This should be compared to a difference of less than 2% for 



these gauges when calculated with the Curtis- Pennington vertex [11]. Even with the CP 



vertex, the variation is significantly greater for covariant gauge choices outside this range, 



particularly for negative gauges [jTq| . Extensions of this work to include nonperturbative 



renormalization were first performed numerically in Refs. [|H|-|T9"||. In these latter works an 
obvious gauge covariance violating term, arising from the use of cut-off regularization and 
present in |TT] , P^[ was omitted (see Refs. [0,0 for a discussion of this). Without the gauge 
covariance violating term the variation of the critical coupling near the Landau gauge was 
again found to be rather small (less than 3% when going from Landau gauge to £ =0.5 [Jill). 

Clearly, the gauge dependence of the (physical) critical coupling is decreased, but not 
eliminated, through the use of a photon-fermion proper vertex which satisfies the WTI. 
In other words, the choice of a vertex satisfying the WTI is a necessary but not sufficient 
condition in order to ensure the full gauge covariance of the Green functions of the theory 
and the gauge invariance of physical observables. The question that arises is whether the 
remaining gauge dependence in the critical coupling is primarily due to limitations of the 
vertex itself or whether it is due to the use of a UV cut-off regulator in these calculations. 

Bashir and Pennington fl"3| , |T4ll have pursued the first of these alternatives and have ob- 
tained, within a cut-off regularized theory, further restrictions on the transverse part of 
the vertex which ensure by construction that the critical coupling indeed becomes strictly 
gauge independent. It is not clear that adjusting the vertex to remove an unwanted gauge- 
dependence is the most appropriate procedure when using a gauge-invariance violating reg- 
ularization scheme. It should be remembered that a cutoff in (Euclidean) momentum, apart 
from breaking Poincare invariance, indeed also breaks gauge invariance. It is precisely 
because of this lack of gauge invariance that UV cut-off regulators have not been used 
in perturbative calculations in gauge theories for many years. Rather, the most common 
perturbative method of regularization in recent times has been that of dimensional regular- 
ization where gauge invariance is explicitly maintained. 
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In this work we report on a DSE study of quenched nonperturbative QED using di- 
mensional regularization in the renormalization procedure. Some early exploratory studies 
have been carried out in the past ||21|| , but to our knowledge this work is the first complete 
nonperturbative demonstration of dimensional regularization and renormalization. Nonper- 
turbative renormalization is performed numerically, in arbitrary covariant gauge, using the 
procedure first developed and applied in Refs. |17|-[19|. In these works the vertex used was 
that of Curtis and Pennington and so, as it is our aim to compare to previous results ob- 
tained with the use of cut-off regularization, we also use this vertex in the current work. 
In quenched QED there is no renormalization of the electron charge and the appropriate 
photon propagator is just the bare one. The resulting nonlinear integral equation for the 
fermion propagator is solved numerically in D = 4 — 2e < 4 Euclidean dimensional space. 
Successive calculations with decreasing e are then extrapolated to e = 0. 

The organization of the paper is as follows: The renormalized and dimensionally regular- 
ized SDE formalism is discussed in Sec. This is followed by some representative numerical 



results in Sec. III. We present conclusions and an outlook in Sec. IVl An appendix details 



the final form of the fermion self-energy equations in D-dimensional Euclidean space. 



II. FORMALISM 

In this section we provide a brief summary of the implementation of nonperturbative 
renormalization within the context of numerical DSE studies. We adopt a notation similar 
to that used in Refs. fT7HT9f, to which the reader is referred for more detail. The formalism 



is presented in Minkowski space and the Wick rotation into Euclidean space can then be 
performed once the equations to be solved have been written down. It is important to 
note that although we use dimensional regularization, we can not make use of the popular 
perturbative renormalization schemes which are usually used in connection with this, such 
as MS or MS. The reason, of course, is that these schemes can only be defined in a purely 
perturbative context. 

The renormalized inverse fermion propagator is defined through 

S-\^p) = A(/i;p 2 ) i- B^- P 2 ) = Z 2 (^e)[^-m Q (e)] -V(ji,e;p) 

= ]f-m(p) -£(//; p) , (1) 

where /i is the chosen renormalization scale, m(/i) is the value of the renormalized mass at 
p 2 = fi 2 , mo(e) is the bare mass and ^(/i, e) is the wavefunction renormalization constant. 
Due to the WTI for the fermion-photon proper vertex, we have for the vertex renormalization 
constant Zi(fi, e) = Z 2 (fi, e). The renormalized and unrenormalized fermion self-energies are 
denoted as p) and e; p) respectively. These can be expressed in terms of Dirac and 
scalar pieces, where for example 

£'G",e;p) = £^(>,e;p 2 ) + X' s (/i,e;p 2 ) , (2) 

and similarly for £(/i;p). We shall for notational brevity not explicitly indicate the depen- 
dence on e of the renormalized quantities A(/i;p 2 ), B(fi;p 2 ) and since for these and 
other renormalized quantities we will always be interested in their e — > limit. The renor- 
malized mass function M(p 2 ) = B(fi;p 2 )/A(n;p 2 ) is renormalization point independent, 
which follows straightforwardly from multiplicative renormalizability [fij 
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The renormalization point boundary condition 



S (p;p) 



— m(/x) 



(3) 



implies that A(fi; /x 2 ) = 1 and m(/i) = M(/i 2 ) and yields the following relations between 
renormalized and unrenormalized self-energies 



S djS (/x; p 2 ) = e; p 2 ) - ^(/x, e; /i 2 ) . 

Also, the wavefunction renormalization is given by 

Z 2 (n,e) = 1 + E' d (/x, e; /x 2 ) 

and the bare mass mo(e) is linked to the renormalized mass m(/i) through 

m (e) = [m(/i) - £' s (/x, e; /?)] /Z 2 (p, e) . 

It also follows from MR that under a renormalization point transformation /i — > fi f , m(p') 
M(n' 2 ) and Z 2 (/x' ', e) = /i 2 ) Z 2 (/x, e) as discussed in Ref. flTj| . 
The unrenormalized self-energy is given by the integral 



(4) 



(5) 



(6) 



r d D k 

V(j*,e;p) = iZi(ji,e)[e(ji)i/\ 2 J ^^^(^ k)T u ([i; k,p)D Xv (p.;p - k) , 



(27T) 



(7) 



where v is an arbitrary mass scale introduced in D dimensions so that the renormalized cou- 
pling e(p) remains dimensionless. Since we are here working in the quenched approximation 
we have Z 3 (p, e) = 1, e = e(/i), and the renormalized photon propagator D^^p; q) is equal 
to the bare photon propagator 



LT{q) 



(8) 



with £ being the covariant gauge parameter. Finally, T u (p; k,p) is the renormalized photon- 
fermion vertex for which we use the CP Ansatz, namely (q = k — p) 



k,p) = rg c (/x; k,p) + T-ete fcV, g 2 ) 7 V - k 2 ) + (p + fc)^ 



(9) 



where Tbc is the usual Ball-Chiu part of the vertex which saturates the Ward-Takahashi 
identity 



K c (^k,p) = -[A(fi;k 2 ) +A(n;p 2 )]Y 



(10) 



(k + pY 
k? — p 2 



[Afa k 2 )-A{^p 2 



+ 



[B(p; k 2 )-B{^p 2 )] 



and the coefficient function t?, is that chosen by Curtis and Pennington, i.e., 



t 6 (ji; k 2 ,p 2 ,q 2 ' 



-~L4(/i; k 2 )-A^;p 2 )}/d(k,p) 



(11) 



where 



fe)= M^ffi±«W. (12) 

The unrenormalized scalar and Dirac self-energies are extracted out of the DSE, Eq. (^), 
by taking |Tr of this equation, multiplied by 1 and p'/p 2 , respectively Note that we use the 



conventions of Muta 22 



7 M 7„ = D YYlf, = (2 - D)Y 

Tr [W] = Tr[l] = 4 

9% = D 

for the Dirac algebra. 

The integrands appearing in Eq. (0) only depend on the magnitude of the internal 
fermion's momentum k 2 as well as the angle 9 between the fermion and photon momen- 
tum. Hence the D-dimensional integrals reduce to 2-dimensional ones, i.e., 

J d D kf(k 2 ,p 2 ,k-p) = J dtt - 1 x J dkk D ~ x J d9 sin D - 2 9 f(k 2 ,p 2 , k ■ p) (13) 



where 



/ 



2ir D / 2 



is the surface area of a D-dimensional sphere. Furthermore, it is possible to express all the 
angular integrals in terms of a single hypergeometric function so that it is only necessary to 
do one integral numerically. The final form of the regularized self-energies is presented in 
the appendix. 

The momentum integration is done numerically on a logarithmic grid and the renor- 
malized fermion DSE solved by iteration. Note that the momentum integration extends to 
infinity, necessitating a change in integration variables. A convenient choice of transforma- 
tion is 



y = y\o {j—j.) > ( 15 ) 

where y\ Q is some lower integration bound and the integration variable t ranges from — 1 to 1. 
The infinite range of the integration also requires an extrapolation of A(fi;p 2 ) and B(n;p 2 ) 
above the highest gridpoint. We check insensitivity to this extrapolation by comparing 
results obtained with a number of different extrapolation prescriptions. In addition, we use 
grids which extend some 20-30 orders of magnitude beyond what is usually used in cut-off 
studies. In summary, we believe that we have verified that the effect of the extrapolation to 
infinity is well-controlled. 
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III. RESULTS 



We present here solutions for the DSE for two values of the coupling a = e^/Aw, namely 
a = 0.6 and a = 1.5. These were chosen so that they correspond to couplings respectively 
well below and above the critical coupling found in previous UV cut-off based studies. The 
gauge parameter is set at £ = 0.25, the renormalization point which we used is /i 2 = 10 8 and 
the renormalized mass is taken to be m(/i) = 400. Note that all results are quoted in terms 
of dimensionless units, i.e., all mass and momentum scales can be simultaneously multiplied 
by any desired mass scale. 

Figs. [IJ and |2| show a family of solutions with the regulator parameter e decreased from 
0.08 to 0.03 for the two values of the coupling. We see that the mass function increases in 
strength in the infrared and tails off faster in the ultraviolet as e is reduced or a is increased. 

Furthermore, it is important to note the strong dependence on e, even though this 
parameter is already rather small. As one would expect, the ultraviolet is most sensitive 
to this regulator, however even in the infrared there is considerable dependence due to 
the intrinsic coupling between these regions by the renormalization procedure. This strong 
dependence on e should be contrasted with the situation in cut-off based studies where it 
was observed that already at rather modest cut-offs (A 2 « 10 10 ) the renormalized functions 
A and M had reached their asymptotic limits. At present it is not possible to decrease e 
significantly below the values shown in Figs. [I] and because of limitations due to numerical 
noise. In order to extract the values of A and M in four dimensions we therefore need to 
extrapolate to e = 0. More sophisticated numerical techniques are being investigated and 
will hopefully allow explicit calculations at smaller e values in the future. 

An extrapolation such as this always involves an added uncertainty in the final result. 
It is fortunate that it is possible to estimate this uncertainty by making use of the fact that 
in the limit e —>■ the renormalized quantities should become independent of the arbitrary 
scale is, which was introduced to keep the coupling a dimensionless in D dimensions. In 
Fig. [3] we show A(p 2 ) and M(p 2 ) evaluated with a = e^/An = 1.5 in the infrared (at p 2 = 1) 
as a function of e for a range of values of v. The results at e = are extracted from cubic 
polynomial fits in e. As may be observed, the agreement between the different curves at 
e = is excellent, being of the order of 0.2%. 

In Fig. |] we show the results extrapolated to e = as a function of the momentum 
(again for a = 1.5). Also shown, although hardly distinguishable, is the result for these 
curves as obtained in the modified UV cutoff based studies, which used a gauge-covariance 
fix to remove an obvious part of the gauge dependence induced by the cutoff |TB[]. Again 
the agreement is very good for a wide range in p 2 and is. Only in the ultraviolet region 
(above say p 2 = 10 12 ) do differences between the curves become discernible. We see in Fig. |] 
that the UV cutoff result is almost indistinguishable from the v = 100 result at essentially 
all momenta. The discrepancies for the v — 1 and 10 cases are greater in the UV due 
to the fact that in this region the errors introduced by the extrapolation procedure in e 
become comparable (for a cubic fit) to the functions' values there. This is supported by the 
following two observations. Firstly, from Fig. [3] we see that the e-dependence is almost linear 
for v — 100, whereas it clearly deviates from linearity for the other cases. Secondly, there is 
little change in the v = 100 extrapolation as we change the order of the fit polynomial up 
to order e 6 , whereas the v — 1 and 10 results show greater variation. We conclude that the 
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v = 100 results in Fig. |3| provide the most reliable e — > extrapolation. It should also be 
noted that the oscillatory behaviour in the mass function first noticed in [17| is reproduced in 
this work, and the results for the unmodified UV cutoff disagree with those from dimensional 
regularization. 

Finally, we present in Tables | and [IT] the bare mass rao(e) and the wavefunction renor- 
malization for the two cases studied as a function of the regularization parameter e and the 
gauge parameter £. It should be noted that these two quantities are by their very nature 
sensitive to the behaviour of A(p 2 ) and B(p 2 ) in the ultraviolet. In particular, while any 
renormalized quantities were found to have a negligible dependence on the precise form of 
the extrapolation of the integrands beyond the highest gridpoint, it was found that Z 2 ({i, e) 
and mo(e) do show some small dependence (< 1%) in our calculations. The values for these 
quantities listed in Tables | and [TI| were obtained by assuming a simple power-behaviour 
of the integrands beyond the highest gridpoint. As was seen in the UV cutoff studies, the 
wavefunction renormalization (for £ > 0) actually decreases as the regularization is removed. 
In the same way the behaviour of Z 2 as a function of the gauge parameter is qualitatively 
the same as observed in the cut-off studies, i.e. at fixed e it decreases as one moves from 
the Landau towards the Feynman gauge. The bare mass mo, on the other hand, appears to 
show a different behaviour to before, at least for large couplings: while in the present work 
it decreases as the gauge parameter is increased for all e, in cut-off studies it only did this 
for moderately small cut-offs. It is quite possible that this different qualitative behaviour 
reflects the fact the even the lowest e's which have been reached here "correspond" to rather 
modest cut-offs. In connection with this note that at the values of e shown in Tables | and [D 
the mass functions in Figs. [I] and|^ are still positive everywhere, the oscillations only setting 
in as one extrapolates toward e — > 0. 
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TABLES 



TABLE I. Renormalization constant /^(/t, e) and bare mass mo(e) as functions of regulariza- 
tion parameter e for a = 0.6 in various gauges £. All solutions are with renormalization point 
/i 2 = 1.00 x 10 8 and renormalized mass m{ji) = 400.0 



e 




m (e) 




£ = o 


£ = 0.25 


£ = 0.5 




e = o 


£ = 0.25 


£ = 0.5 


0.08 


1.000 


0.957 


0.917 




2.30xl0 2 


2.29xl0 2 


2.28xl0 2 


0.07 


1.000 


0.944 


0.891 




1.91xl0 2 


1.90xl0 2 


1.89xl0 2 


0.06 


1.000 


0.924 


0.854 




1.46xl0 2 


1.45xl0 2 


1.44xl0 2 


0.05 


1.000 


0.896 


0.802 




9.75X10 1 


9.66X10 1 


9.57X10 1 


0.04 


1.000 


0.851 


0.725 




5.07X10 1 


5.00X10 1 


4.94X10 1 


0.03 


1.000 


0.779 


0.606 




1.59X10 1 


1.56X10 1 


1.53X10 1 


TABLE II. As above, but for a = 1.5 


e 




m (e) 




e = o 


f = 0.25 


£ = 0.5 




£ = o 


C = 0.25 


£ = 0.5 


0.08 


1.000 


0.897 


0.805 




9.20X10 1 


9.06X10 1 


8.92X10 1 


0.07 


1.000 


0.865 


0.749 




5.51X10 1 


5.40X10 1 


5.28X10 1 


0.06 


1.000 


0.821 


0.674 




2.58X10 1 


2.51X10 1 


2.44X10 1 


0.05 


1.000 


0.759 


0.576 




7.94x10° 


7.60x10° 


7.28x10° 


0.04 


1.000 


0.669 


0.447 




1.10x10° 


1.02x10° 


9.54X10" 1 


0.03 


1.000 


0.535 


0.286 




2.63xl0~ 2 


2.27xl0" 2 


1.96xl0" 2 
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IV. CONCLUSIONS AND OUTLOOK 



We have reported here the first detailed study of the numerical renormalization of the 
fermion Dyson-Schwinger equation of QED through the use of a dimensional regulator rather 
than a gauge invariance-violating UV cut-off. The initial results presented here are encour- 
aging. Firstly, we have explicitly demonstrated that the approach works and is independent 
of the intermediate dimensional regularization scale {u) as expected. Secondly, we have 
seen the interesting result that our calculations using dimensional regularization agree with 
modified UV cut-off calculations within the current numerical precision, but disagree with 
the unmodified UV cutoff ones. 

A significant practical difference between the dimensionally and UV cut-off regularized 
approaches is that in the former it is at present necessary to perform an explicit extrapo- 
lation to e = whereas in the latter it was found that for a sufficiently large choice of UV 
cut-off the results became independent of the cut-off. This need to extrapolate, together 
with the high precision that one needs to attain in order to make meaningful comparisons 
with cut-off based studies, makes numerical dimensional regularization and renormaliza- 
tion of Schwinger-Dyson equations a rather formidable task. Nevertheless, although we 
are presently investigating whether it is numerically possible to extend the studies to even 
smaller values of e in order to improve the precision of the e — > extrapolation, the extrapo- 
lation to e = appears to be well under control, at least for values of the fermion momentum 
away from the ultraviolet region. 

Having demonstrated the numerical procedure of renormalization using dimensional reg- 
ularization, we now plan to study chiral symmetry breaking and in particular hope to extract 
the critical coupling as a function of the gauge parameter. Results of this ongoing work will 
be presented elsewhere. The eventual aim is to extend this treatment to the case of un- 
quenched QED and to a systematic study of electron-photon proper vertices. 
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APPENDIX A: FINAL FORM FOR THE REGULARIZED FERMION 

SELF-ENERGIES 

In the quenched approximation, all angular integrals in the Dirac and scalar regularized 
self-energies defined by the Euclidean analogue of Eq. (|7]) may be expressed in terms of the 
integrals 

4» = - L^f) - Fd9 (^±y) n sin D ^e for n = -1,1, 2 (Al) 
1 ( — y~ ) V 71 " J0 \ z / 

where we have defined dimensionless quantities x = p 2 /v 2 , y = k 2 /V 2 , z = q 2 ' jv 2 = x + y — 
l^fxy cos# and w—y/x. Similarily, we shall for convenience define dimensionless versions of 
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A(fi 2 ;p 2 ) and B(fi 2 ;p 2 ), namely a(x) = A(fi 2 ;p 2 ) and ub(x) = B(fx 2 ;p 2 ) (we suppress the 
dependence on fx here in order to make the notation less cumbersome). Explicit evaluation 
(for < w < 1) yields 



I?(w) = (l + w) 2 F 1 (l,e;2-e;w) 



(A2) 



h D H = 2 ( T 

and for w > 1 one may use 



1+w 



w 



1 - 2e 



jfH + 1 - e) , 



Defining 



d E (x,y) 
Aa(x,y) 



d E (x,y) 



x + y 
a(x) - a(y) 
(1 -w)/(l + w) 



(A3) 

(A4) 
(A5) 



(and similarly for Ab(x,y)) the regularized self-energies in the quenched approximation 
become 



S d (^,e;x) = — (4tt) 



2/~ 



x 



4vr v y r(2-e)7o V 2 (?/) + &*(?/) 

^)^(i- £ )(i-/fH) 
(| + (i-0(i-^)- 



+a(y) 



2 

Aa(x,y) 



3-2e /1-iu 



2d E (x,y) \l + w 



+-^-AS(x,y) 



(1 + w) 2 

(5 + (! - (i - c)) (A» - 1) + 



(A6) 



and 



Z s (ix,e;x) = — (4tt) 



,1-c 



4vr y r(2-e)7 ya 2 (y) + 6 2 (y)a;(l +w) 
x | 6fe) aW+aM (3 _ 2f + 0/i>) 



-a(y)Ab(x,y) 
+b(y)Aa(x,y) 



(l + i 1 - (i - 0) (A» - 1) + if^^H 



£ 3-2e /l-u^ 2N 



2 ' 2^(x,y) U + w 



/f( W ) 



(A7) 
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It should be remembered that the equivalent expressions in cut-off regularized quenched 
QED suffer from an ambiguity due to the lack of gauge invariance in those calculations: one 
obtains different results at this stage depending one whether or not one has made use of the 
Ward-Takahashi identity in the initial stages of the calculation |T^Jl3| . The difference shows 
up in the terms multiplied by a 2 (y) and b 2 (y). Readers may readily convince themselves 
that no such ambiguity exists in the present work. In order to do this, the following identity 
is of use: 







Jo \ y-x 



2e- 1 



y-x 



x + y 



1 




) 



(A8) 
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FIGURES 

oc=0.6 solutions for different £ 
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FIG. 1. The finite renormalization A(p 2 ) and mass function M(p 2 ) for various choices of 
the regulator parameter e. These results have coupling a = 0.6, gauge parameter £ = 0.25, 
renormalization point fi 2 = 10 8 , renormalized mass m(n) = 400 and scale v = 1. In the low p 2 
region the smallest e has the largest value of M{p 2 ). 
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oc=1.5 solutions for different £ 
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FIG. 2. Same as Fig. 1 for a = 1.5 
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oc=1.5 extrapolation to £=0 at p =1 
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FIG. 3. The finite renormalization A(p 2 ) and mass function M(p 2 ) evaluated at p 2 = 1 for 
various values of the regulator parameter e and extrapolated to e = by fitting a polynomial cubic 
in e. Shown are results for three scales v = 1, 10 and 100. All other parameters are those of Fig. 
2. The different scales coincide at e = to an accuracy of approximately 0.2%. 
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a=1.5 solutions extrapolated to 8=0 
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FIG. 4. The finite renormalization A(p 2 ) and mass function M(p 2 ) extrapolated to e = at 
every momentum point for three different scales v = 1, 10 and 100. They were calculated by fitting 
a cubic polynomial in e at each momentum point. All other parameters are those of Fig. 2. Also 
shown is the result obtained using a modified UV cutoff. The UV cutoff and extrapolated v = 100 
results are indistinguishable in these plots. 
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